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We compute the nonperturbative decay of supersymmetric flat directions due to their D-term 
potential. Flat directions can develop large vacuum expectation values (vevs) during inflation, and, 
if they are long-lived, this can strongly affect the reheating and thermalization stages after the 
inflation. We study a generic system of two 17(1) or SU(2) flat directions which are cosmologically 
evolving after inflation. After proper gauge fixing, we show that the excitations of the fields around 
this background can undergo exponential amplification, at the expense of the energy density of the 
flat directions. We compute this effect for several values of the masses and the initial vevs of the two 
flat directions, through a combination of analytical methods and extensive numerical simulations. 
O |' For a wide range of parameters the flat directions decay within their first few rotations. 

PACS numbers: 98.70.Cq 

I. INTRODUCTION 
t-H ■ 

Flat directions are prevalent in the MSSM [l[ and in its simplest extensions. Scalar fields can develop a large vev 
' along a flat direction during inflation [2|,|3j], though in supergravity models, this typically requires a nonminimal Kahler 
fSJ . potential 0, H, Q. If large field values are excited, flat directions can have several important effects in cosmology Q. 
' For example, in the presence of a non- vanishing baryon number violating operator, the flat direction may be associated 
, with a finite baryon number density and could account for the production of the baryon asymmetry of the Universe 
through the Affleck-Dine mechanism [H, Q . More recentl y, it has been argued that the flat directions may also affect 
OO ! the thermal history of the universe after inflaton decays There are two main effects. First there is the delay of 
■ the thermalization of the inflaton decay products (resulting in a low reheat temperature) due to the fact that the vev 
^ \ of the flat directions provides a large mass to gauge bosons, which in turn suppresses the rates of the processes needed 
• , for thermalization [lCj] • 1 The second effect is that the oscillations of the flat directions can come to dominate over 
the (relativistic) inflaton decay products @ , so that the relevant stage of reheating is associated with the decay of the 
^ , flat directions, rather than that of the inflaton field. 

For definiteness we assume that there is only a single mass scale associated with the inflaton potential, set by the 
COBE normalization [12] to be about ~ 10~ 7 M p (where M p is the Planck scale). We also assume that the 
inflaton couples to the visible sector only through gravity. It was shown in p"l"j that, under these circumstances, the 
delay of the thermalization of the inflaton decay products takes place only if the vev of the flat direction generated 
during inflation satisfies (f>o ^ ev 3 / 2 Mp^rn^/m^p 2 , where a is a characteristic coupling strength between fermions and 
the gauge bosons (hence, of order of the gauge fine structure constant) and the mass of the flat direction. For 
a 2 ~ 10 -3 , and for ~ 100 GeV (since the mass comes from soft SUSY breaking terms) this gives 4>o ~ 10 17 GeV . 
It was also shown that the energy density of the flat directions comes to dominate over that of the inflaton decay 
products for </> ~ Mp^ 3 m 5 J 12 /m 3 ^ 4 ~ 10 17 GeV. 

For our canonical choices of masses, the two conditions turn out to be the same. That is, the vevs of the flat 
directions need to be very large to produce a delay in thermalization or for the flat directions to dominate the energy 
density. A limit on the vev of flat directions may be derived from higher-order superpotential terms of the form 
4> n /M n ~ 3 , with n > 3 . Gauge invariance can protect MSSM flat directions and for each class of flat directions 
there is a minimum value of n for which the flat direction is lifted No MSSM flat direction is protected beyond 



1 As pointed out in , this delay argument assumes that the vevs of the flat directions break all the gauge symmetries of the Standard 
Model. If this is not the case, some of the gauge bosons remain light, and thermalization may proceed through them. 
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n = 9 . Clearly, the fact that such terms are allowed by gauge invariance, does not imply that they are indeed in the 
superpotential, since their presence is model-dependent (they may for instance be forbidden by additional symmetries) . 
These non-renormalizable terms, if present, can be used to determine the vev along a flat direction during inflation. 

Given sufficient time, a vev of order </>o ~ (m^ Af ™~ 3 ) 1 ^" ^ is generated, though typical inflationary models built on 
a single scale will produce a vev of order O 2 ~ H 3 t ~ m^M P Q. In order to obtain a vev as large as > 10 17 GeV, 
we must require that nonrenormalizable terms with low n are absent. For M — M p , and = 100 GeV, we find 
that all terms with n < 11 must be absent. It should be noted, however, that loop corrections in supergravity during 
inflation may lead to much larger vevs, particularly in models of the no-scale form 0, @. 

Delayed thermalization and/or domination by flat direction requires that the flat directions themselves are long- 
lived. As pointed out in this appears to be the case if one only considers the perturbative decay of the flat 
directions. The perturbative decay rate of flat directions is suppressed by their large vev, ~ m^/cj) 2 Q. For 
0o ^ 10 17 GeV , one finds indeed that the perturbative decay of the flat directions takes place after they dominate. 
More specifically, one finds ~ m 5 / 3 /Mp^ 3 [11]. Since the typical timescale for an oscillation of a scalar field along 



a flat direction prior to decay is m^ 1 , one finds that flat directions perform about (M p /m^,) 2 ' 3 ~ 10 10 oscillations 
before their perturbative decay. 

However, it may be possible that the decay of a flat direction vev is controlled by nonperturbative effects. This 
possibility was first discussed in [l3|, fl4| . where it was assumed that the scalar field along a flat direction has an 
interaction of the type AV = <? 2 |0| 2 |x| 2 with another scalar field x 2 ■ The cosmological evolution of the flat direction 
then results in a time dependent effective mass for x- I n principle, this interaction could lead to a nonadiabatic 
variation of the frequency for x, namely 6j x /uj x > 1, and, consequently, to the nonperturbative production of x &t the 
expense of the energy stored along the flat direction. Such an effect is often included for the inflaton field, leading to 
nonperturbative inflaton decay or preheating [16| . It was, however, found that this effect is absent in the case of flat 
directions. The reason being, that contrary to the case of inflation, flat directions are complex fields that typically do 
not pass through \<fi\ = during their "oscillations". 3 Therefore, the frequency cu x is always large enough to ensure 
that u) x jio x < 1 at all times. 

Although formally correct, this conclusion heavily relies on the assumed coupling AV^ = g 2 |0| 2 |x| 2 - However, 
the couplings of MSSM flat directions to other fields is much more complicated than this, and may give rise to 
nonperturbative effects not considered in [l3l. In particular, the D-term potential results in a large nondiagonal 
and time dependent mass matrix for several MSSM fields which depends on the (time dependent) background of 
the flat directions. Even if the eigenvalues of such a mass matrix are constant, or slowly evolving, this mass matrix 
"rotates" (in field space) due to the rotation of the flat directions in their complex plane. Such an effect can result 
in a very strong (exponential) production of MSSM fields, at the expense of the energy density stored in the flat 
directions, which therefore would decay on a much quicker timescale (O (10 — 100) rotations) than that suggested by 
the perturbative decay rate (Tlj . 

It was shown in [ll| that the non-perturbative decay of a flat direction vev requires the presence of two or more 
simultaneous flat directions. In the case of a single direction, the rotation of the flat direction in its complex plane can 
be absorbed by a gauge redefinition. This should not be seen as a serious problem, since it is easy to verify that only in 
exceptional cases does a single flat direction exclude all others from being excited. Therefore, if the conditions in the 
early universe allowed the formation of a large vev for a flat direction (sufficiently long phase of inflation, a nonminimal 
Kahler potential, and absence of the nonrenormalizable interactions of low orders) , then one should conclude that a 
system of non-mutually exclusive flat directions developed large vevs. For the case of multiple flat directions, it was 
shown [ll| that for several simple cases the number of degrees of freedom involved is sufficiently large that the mixing 
in field space remains even after all the gauge redundancy has been removed. However, the actual computation of 
the decay was performed only in toy models, characterized by a potential for the fields that has the same structure 
of the MSSM D-term potential: in such models, the flat directions experienced very rapid nonperturbative decay. 
The computation of the decay in complete models, with all the gauge interactions properly taken into account was 
postponed to a separate publication, and that is the main goal of the present work. As we will see, the result of the 
detailed calculation presented here agrees with those of (ill ]. 

The paper is organized as follows. In Section [III we introduce U(l) and SU(2) models with one and two flat 



2 Nonperturbative effects in the presence of supersymmetric Q-balls were discussed in [T3| . 

3 This is due to terms which break the radial symmetry of the potential of the flat direction, and which therefore generate an angular 
motion. One can however envisage a situation in which such terms are absent, so that the motion is purely radial. In this case, preheating 
effects like parametric resonance [TJ| or instant preheating [l7ll lead to the fast decay of the flat direction [18| . 
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directions, and we motivate why multiple flat direction models are necessary in order to undergo nonperturbative 
decay. In Section HTT1 we outline the formalism for computing this decay, discuss under what conditions a large effect 
can be expected. In Section IIVI we perform a detailed analysis of the models introduced in Section [TTJ This study 
allows us to clarify why the nonperturbative production that we are studying is absent in the single flat direction 
case, but is instead present for multiple flat directions. In Section [V] we present our numerical results for the particle 
production when two flat directions are excited. The computation is summarized in the discussion Section PVTl where 
we also compare our results with others that have appeared in the literature after The paper is concluded by a 
number of Appendices, in which we collect some of the technical steps of our calculations. 



II. MODELING NONPERTURBATIVE FLAT DIRECTION DECAY 



In this section, we present the models considered in this work. These are specific models of gauged U(l) and SU(2) 
flat directions. We derive the light degrees of freedom in the unitary gauge which can play a role in the non-perturbative 
decay of the fields excited along the flat directions. 

We first consider a U(l) gauge theory with two fields of charge +e and — e, with Lagrangian 4 

£=\Dfa\ 2 + \D ( f> 2 \ 2 -jF 2 -V(4> 1 , fa) 



Df, = ( <9 M - — qiAp ) <j)i (1) 



where qi 2 = 1, —1 . The dominant potential term, originating from a D-term, is 

V Dl = j(|<£i| 2 -|<fc| 2 ) 2 (2) 

which admits the flat direction \<j>i\ = . In the unitary gauge, this theory has 6 degrees of freedom: the three 
polarizations of the gauge boson, two light modes, corresponding to the real and imaginary part of the excitations 
along the flat direction, and one massive Higgs. 

This decomposition is worked out in detail in appendix [AJ We expand the scalar fields about some background 
value <&, for which the D-term potential @ vanishes: <\>\ — $ + 8<f>x and <fi 2 = $ + 8<p 2 - Inserting these expansions in 
the potential one finds that the light modes are defined by 8<pi + 8<p2- These fields can be rotated into the real and 
the imaginary excitations of the flat direction itself. They are decoupled from the other fluctuations at the linearized 
level, and their eigenfrequencies vary only adiabatically with time. Therefore, such fields do no undergo parametric 
resonance. 

The transverse components of the vector fields are also decoupled at the linearized level. On the contrary, the 
longitudinal component is coupled to the Higgs field. However, as we show below, this coupling is strongly suppressed 
by the smallness of the flat direction mass (with respect to its amplitude) and can be neglected in the physically 
relevant cases. Therefore, all four of the massive fluctuations can be thought as decoupled from each other. Their 
mass terms are of the form 

VD _I e 2 |$| 2 |^| 2 + ie 2 |$|%^ (3) 

where the Higgs fluctuation Su is obtained from the combination 8<j>i — 8<j>2 (see below). The equality of the mass 
of the Higgs and of the gauge boson is due to the fact that the model is chosen to reproduce the bosonic sector of a 
supersymmetric theory. Actually, the masses of these fields receive additional contributions from the soft supersym- 
metry breaking mass terms m 2 \4>i\ 2 of the two scalars, which are however much smaller than the amplitude |$| and 
can be ignored for now. These terms break the degeneracy between the masses of the gauge boson and the Higgs. 
The potential ((3]) leads to the time dependent effective mass of 8h and 

m 2 cS = e 2 \^(t)\ 2 (4) 

In analogy to what has been computed for inflationary preheating [l6[ , one can investigate if the time variation of the 
mass is quick enough to lead to nonperturbative particle production. This nonperturbative decay of the background 



4 The normalization for the electric charge is chosen to match that of the SU(N) generators. 
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fields takes place whenever the frequency of the quanta varies nonadiabatically, u> > uj 2 . For relatively small momenta 



p, one finds uj 2 — p 2 + m 2 s ~ e 2 1 <I> (t) \ 2 , and thus 



— ~ — (5) 
uj 2 |$| (b) 

using |$| ~ m|$|, since this is the scale that governs the dynamics of the background. Since, m is of order the 
electroweak scale, while $ is close to the Planck scale.it is difficult to construct a situation in which this quantity 
becomes near unity. Based on this consideration [ID, [l4[ , it was concluded that the flat direction decays only 
perturbatively, leading to a suppressed reheating temperature and a solution to the gravitino problem. However, this 
conclusion strongly depends on the coupling © assumed in these analyses, and in the resulting mass term In 
several concrete cases, the coupling is more complicated. 

This is indeed the case if more flat directions are present. To see this, we add two more fields to the Lagrangian 
also of charge +e and — e (other choices are possible as well), and write the D-term contribution to the potential, 



v D2 = ^(I0i| 2 -|0 2 | 2 + |0 3 | 2 -|04| 2 ) 2 (6) 

We discuss the case in which two flat directions are present in this model (notice that three flat directions are also 
possible). There are ten real modes around these directions (coming from the four complex scalars and the gauge 
boson); in the unitary gauge, they can be decomposed as follows (see appendix [51 for the detailed computation): three 
are associated with the massive vector, four provide shifts to the two flat direction vevs, and the three remaining 
degrees of freedom (one heavy and two light ones) are coupled to each other in a Higgs mass matrix. The longitudinal 
component of the vector field is also coupled to these fields, so this is actually a system of four coupled fields; however 
- analogously to what happens for the single flat direction case - this coupling turns out to be negligible in the regime 
of small masses vs. amplitudes that we are interested in. 

Denoting the flat direction vevs by $ and $, the Higgs mass matrix possesses a heavy part proportional to e 2 $ 2 
(and e 2 $ 2 ) which originates from the D-term ([5]), and a light part composed of soft supersymmetry breaking masses 

m|. as well as terms of order |4>/$| 2 and |$/$| 2 which appear from gauge fixing and from making the kinetic terms 

canonical. The light terms also come with time dependent prefactors such as , l£l _ and , ffl = . The resulting 

Vl*l 2 +I*l 2 Vl*l 2 +I*l 2 

time dependence of the Higgs mass matrix in general cannot be removed by a field redefinition, and so all of the 
eigenmasses retain this time dependence. The spectrum consists of one heavy state and two light states of order m t j >i 
and in the limit of zero momentum, they all oscillate at the frequency . Thus, for momenta lower or comparable 
to the supersymmetry breaking masses, 

UUgM _ j (?) 
W light 

which suggests a nonadiabatic evolution of these light states. Similarly, it can be shown that the corresponding 
eigenvectors of these light states Cu g ht will execute relatively fast rotation in the space they span, 

\Clight | ^ ^ ,gs 
flight 

which can also lead to nonadiabatic evolution. The second effect was highlighted in [llj], and it is discussed in detail 
in Sections IIIII and I VII below. Both effects must be considered when there are multiple fields involved in the Higgs 
mass matrix. 

This argument for nonadiabatic evolution relies upon the existence of light eigenstates in the Higgs mass matrix 
and also on their time dependence. The existence of these light states is not uncommon in the presence of multiple 
flat directions. For example, in 11] it was shown that the simultaneous presence of the two flat directions LLE C and 
QLD C results in a Higgs mass matrix with time dependent light states mixed with one heavy state. It is expected 
that the results ([7]) and © will be valid also for this case. 

It must be stressed that the rapid change of the light states is not guaranteed, and in particular it will be suppressed 
when the ratio |$|/|$| is very large or very small. In these limits, there is a suppression in the Higgs mass matrix 

via the above mentioned factors , ^ . and -=J2L= . A suppression of this kind was also noted in [2(1, but a 

Vl*l 2 +I*l 2 ^/|*| 2 +|*| 2 1 — 1 

concrete calculation of the effects was not performed. Our numerical studies indicate that even if the ratio of maximum 
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values of the vevs is as large as ~ 10 +3 , the ratio |$|/|$| is still of order one during some short periods of time. We 
find that even in cases such as this, there is sufficient nonadiabatic evolution during these short periods of time to 
result in the rapid decay of the flat direction. This effect will be demonstrated in Section fVl 

In the remainder of this section, we discuss the background evolution as well as the symmetries of the background 
under changes in the initial conditions and parameters. For the two field case, we consider the potential 

m 2 (\fa\ 2 + I'M 2 ) + K^l +h.c.) + ^(I0i| 2 - I02| 2 ) 2 (9) 

where the last term is the D-term potential ((2]) , while m 2 is a soft supersymmetry breaking mass of order the electroweak 
scale. This potential admits the D-flat direction with vev assignments (fa) = (fa) = <f>. 5 The quartic term 
proportional to A may be present due to one loop contributions to the effective potential and are assumed to have 
magnitude A ~ r|j^p in analogy to the baryon number violating operators in the Affleck-Dine mechanism Q. Here, 
$o is the flat direction vev at the time the flat direction begins its oscillation. These terms provide an initial angular 
motion in the complex plane and become subdominant to the other terms in the action as the universe expands. 
The vev subsequently spirals down slowly to the origin as discussed in detail below. Without the quartic terms, the 
dynamics of the flat direction vev would be a straight line trajectory in the complex plane through the origin. 

The cosmological context for the evolution of the flat directions, is a Universe initially dominated by the scalar field 
oscillations of the inflaton before it decays. Thus initially we assume a matter dominated Universe with the dominant 
contribution being in the inflaton energy density p^. We then have the background equations of motion 



$ + 3f $ + m 2 $ + 2A$* 3 = 



R 3MS 



2|$| 2 + 2m 2 |$| 2 + A($ 4 + h.c.) + (10) 



where we have assumed a Friedmann-Robertson- Walker metric ds 2 = dt 2 — R 2 dx 2 in which R is the scale factor. As 
we will see, it will be convenient to express the flat direction in terms of two scaled real scalar fields, 

* = ^e* (11) 

The equations of motion may be rewritten using conformal time 77 with dt = Rdrj as, 

F" + (m 2 R 2 - SjL - £' 2 ) F + I F 3 cos (4 E) = 
(F 2 F 4 sin (4S)=0 

^(iQ 5 '' = *0-&Vm+Pl>) (12) 

where a prime denotes a derivative with respect to the conformal time, and where the Friedmann equations were 
combined in a way so as to eliminate the kinetic terms for the scalar fields. This form (fT2")) of the background 
equations is used for our numerical analyses shown below. 

The background equations of motion possess two symmetries. Given a solution ^(i) to the equations (TIT)]) with 

initial value $a(0) = $o, and recalling that A cx we see that 7$^(t), with 7 a real constant, is also a solution 

to the equations of motion. Thus, changing the initial value for the field, only changes the solution through a scale 
change. The background equations have a second symmetry under a change of the soft mass m — > \x m, with \i real. 
Now, $ a (fJ- 1) is a solution to the new background equations of motion. Thus changing the mass parameter in the 
equations only changes the timescale on which the solution evolves. Both symmetries will be shown in Section IVT1 to 
translate into approximate symmetries for the decay of the flat directions. 6 



Note, however, that taking mi 5^ m,2 will result in flat direction vevs with more complicated time evolution. 

For this argument, we neglect the subdominant contribution to the Friedmann equation dlOI I from the energy density of the fiat direction. 
The initial energy of the inflaton field is chosen so that the flat direction is initially at rest due to Hubble friction. The flat direction 

1/2 n 

starts evolving when the Hubble parameter oc p ! drops below its mass. Starting with a higher value for simply amounts in waiting 
for a longer time before this evolution starts. Therefore, under the rescaling m —> pm, we also rescale p^ — » p? p^, . 
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This model can be extended to the four fields model with potential 

m 2 (\fa\ 2 + \fa\ 2 ) + rn 2 (\fa\ 2 + \fa\ 2 ) + A(^ 2 + h.c.) + \(<f>jfa 4 + h.c.) 

+ j(\fa\ 2 -\fa\ 2 + \fa\ 2 -\fa\ 2 ) 2 (13) 

characterized by the D-term potential ([S]), soft supersymmetry breaking terms, and quartic interactions. This model 
admits the two flat directions characterized by the vevs ((pi) — (fa) — $ , (fa) — (fa) — $ ■ The background equations 
of motion are an immediate extension of (|12|) ; we simply have two copies of the first two equations in (|12|) (one copy for 
the amplitude and the phase of and one for the amplitude and the phase of $), while the last equation in (fT2|) only 

2 ~ ~ 2 

changes in the sense that the potential for both fields is present on the right hand side. Since A oc jj^p-j and A oc j^-p > 

the system possesses the same symmetries of the previous case under the simultaneous rescalings <£> ^ 7 $, <£> ^ 7 $ 
orrfn/im, m-t/im. 

In this work, we study the decay of the background flat direction(s) into the fields involved in the potentials (0 
and (TlU)) . We do so by studying the linearized theory for the perturbations, and by studying whether the system 
of fluctuations undergoes nonperturbative production due to the coherent motions of the flat direction(s). 7 This 
linearized computation can be extended to nonabclian flat directions of physical interest, as for instance SU(N) flat 
directions. The reason for this is that the differences due to the nonabclian structure only arise at higher than quadratic 
order in the gauge fields. Such interactions are not present at the linearized study of fluctuations, as long as the gauge 
fields involved do not have any background expectation value. 

For simplicity, we focus on the SU(2) case. Consider the model with potential 

m 2 (\fa\ 2 + \fa\ 2 ) + A [(e^faifaif + h.c] + y E (#^1 + ^faf (14) 

a 

where i,j are SU(2) doublet indices, and the r a are the gauge generators T a — \ (where a a are the Pauli matrices). 
This potential admits the flat direction fa = (<&, 0) , fa = (0, $) . As we show in Appendix [C] the physical modes 
around this background consist of the real and imaginary excitations of the flat direction, of three massive gauge 
bosons, and of three Higgses. The quadratic actions for the gauge bosons and the Higgses are three copies of those 
obtained in the U(l) case. 

We can extend this model to contain four complex doublets, with potential 



m 2 (|^| 2 + \fa\ 2 ) + rh 2 (\fa\ 2 + \fa\ 2 ) + A [(e tJ fa t fa b ) 2 + h.c] + A [(e^fa.fa,) 2 + h. 



+ ^Yl + ^r a fa + far a fa + far a fa) 2 (15) 



We can now have two flat directions, with the vev assignments fa — ($, 0) , fa = (0, , fa = ( 0^ , fa — ^0, $j . 

Also in this case, beside the real and imaginary parts of the two flat directions, the physical modes around this 
background consist of three copies of those obtained in the U(l) case; namely, three massive gauge bosons, three 
Higgses, and six light fields. Each Higgs is coupled to one longitudinal boson and two light fields exactly as for two 
U(l) flat directions. 

Due to this correspondence between the linearized U(l) and SU(2) (which can be generalized to SU(N)) systems, it 
is sufficient to explicitly perform the numerical computations only for the U(l) cases. 



III. FORMALISM FOR THE QUANTUM EVOLUTION OF COUPLED SCALAR FIELDS 

Here we outline the formalism of which is necessary to quantify the rate of the nonperturbative decay of flat 
directions. We begin with a Lagrangian with an arbitrary number of real scalar fields 'J' = {ipi,ip2, ■■•'4 , n} and with a 
time dependent mass matrix M 2 (rj), 

C = ^d^ T d^ 1 ^ - i* T M 2 (?7)^ 



7 We comment on the effect of the nonlinear interactions in the Discussion section I VII 



7 



where we are using the conformal time ?y, and derivatives with respect to conformal time are again denoted with a 
prime. One may diagonalize M 2 with a time dependent rotation matrix C{rj), 

C T (ri)M 2 (r])C(r 1 ) ^ m 2 d {r,) diagonal (16) 

The eigenfrequencies of the system are the elements of the diagonal matrix 



cu = \Jk 2 l +m 2 (17) 
where k is the momentum. The fields in the diagonal basis are 

§ = C T * 

so if the mass matrix is evolving slowly enough, the column vectors of C denoted by Ci will be the physical eigenstates 
of the system. The equations for the quantum evolution of the state of this system were shown in [19( to be 

a' = (-iuj-I)a+^-j)(3 

F = {iu-I)0+(£-j y )a (18) 
where a and (3 are matrices of Bogolyubov coefficients for the set of scalar fields, $ and 

i, J = i (v^r-^= ± -^iW) , r = c T C (19) 



Note that T and / are antisymmetric matrices while J is symmetric. In addition, the Bogolyubov coefficients obey 
the constraints, 

aa* - (3*(3 T = 1 (20) 

a/?t - /3*a T = (21) 

and the occupation number of the i'th state can be shown to be, 

ni{rj) = (p* /3 T ).. , no summation on i (22) 

Note that the Gamma matrix may be written in terms of the eigenstates Ci and its elements contain the information 
for the rate of change of the eigenstates, 

= Ci ■ Cj 

C'i = —TjjCj 
j 

There are two qualitatively different parts to the evolution equations (| 18[) ; a "unitary" part represented by the anti- 

Hermitian matrices (±iuj — I) and a "nonunitary" part represented by the symmetric matrix — Whenever 

the nonunitary part vanishes, the total occupation number given by Tr[/3*/? T ] is conserved. However, the matrix 
(3* (3 T itself is not constant in general indicating a conversion from one species of particles to another as a function of 
time. The non-unitary part will boost (or contract) both a.cv and /3*/3 T in (|20p while keeping the difference invariant, 
so this term will change the total occupation number of the system. It will be convenient to rewrite the evolution 
equations (fT5|) by extracting factors of y/u, 

i )=I-j= ) V^a + \[u [ ^2 l -=J-^= ) \/uj(3 



The equations are now written entirely in terms of the following matrices, 



/UJ \/UJ I oj 2 Jlu Juj lu 2 \ lu lu 
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One recognizes the diagonal part of A as the adiabatic parameter from the single-field analysis. The off-diagonal parts 
of this matrix may also be interpreted as adiabatic parameters as follows. One expects nonadiabatic evolution when 
either the oscillation frequency of the physical state or the field composition of the state (its eigenvector) has changed 
faster than its period l/u). The latter condition on the field composition is clearly relevant only if there are two or 
more fields from which the state is composed. Non-adiabatic evolution occurs when the condition (jSJ) is satisfied, 
\C[\ > LUi or more specifically when, 



\ 



r ^ 



This condition is not quite correct however, since in the case of degenerate eigenstates, the eigenvectors have no 
preferred directions in their subspace, and so the rate of change of these eigenvectors in this subspace is not physically 
meaningful. This feature is taken into account by the A matrix (as well as the J matrix) which has the form, 



An = — , no summation 



Aij = Tij ( ) , i 7^ j and no summation (26) 



and it handles both limiting cases, 



lim Aij =0 , i ^ j 
lim Aij = ( — ) , i ^ j 

i+>j>>uii V U!i I 



so for nondegenerate frequencies the condition ([23)) is asymptotically true. We thus take the elements of A as our 
adiabatic parameters. If any one element of the matrix is greater than unity \Aij\ > 1, we expect nonadiabatic 
evolution. We note finally that large A^ are a necessary condition for net production of quanta, but not a sufficient 
condition. 



IV. NONPERTURBATIVE DECAY OF SINGLE OR MULTIPLE FLAT DIRECTIONS 

In this section, we start the study of the nonperturbative decay of one flat direction with potential © and of 
two flat directions with potential ([13]) . The spectra of fields for the two cases are worked out in appendices lAl and 
iBl respectively. The results of these computations are summarized in the two Subsections below. The study of the 
quadratic action for these modes allows one to understand whether the nonperturbative decay of the flat directions (s) 
takes place. The actual numerical computation is presented in the next section. 



A. Single flat direction 

The single flat direction case is characterized by the potential ([9]), and the background values 

= <<fe) = 7^e lS , (A„)=0 (27) 

The perturbations of these fields can be most usefully written as 

1 e iE 
(fyi+<fya = —(S r +iSi) , 6<p! - 84> 2 = -5- (S H + i S G ) (28) 

The mode 5q is the only one in this decomposition that varies under an infinitesimal gauge transformation (since 
the two fields have opposite U(l) charge), and it can be set to zero in the unitary gauge. 
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The two fields 8 r and Si are the physical excitations along the real and imaginary directions of the flat directions. 
They are decoupled from the other fluctuations at the linearized level. However, they are coupled to each other in 
their mass term, due to the quartic term in ([5]), 



J_ f m 2 R 2 - + M F 2 cos (2S) -^F 2 sin(2£) 
lflat ~ B? V -^F 2 sin (2S) m 2 R 2 - 3£ - ^ F 2 cos (2E) 



(29) 



corresponding to the eigenmasses 



2 2 i?" 3 A F 2 

m 1;2 = m "^3 ±^-^2 ( 3 °) 



The time variation of these eigenmasses is adiabatic. Moreover, the off-diagonal quartic term, already subdominant at 
the start, quickly becomes negligible as the flat direction starts evolving. Consequently, the nonperturbative production 
of these fields is negligible as well. 

The remaining mode Sh is the Higgs field of the model. It is coupled to the longitudinal component of the gauge 
boson. The quadratic action for these two modes (after integrating out the non-dynamical fluctuation Aa) is given 
in (|A15|) . We study this action in the phenomenologically relevant case in which the masses of the flat directions are 
much smaller than their amplitudes. More accurately, we work in the limit of 

mR, — (31) 

For bookkeeping, in the following we denote by e the ratio between any of the terms on the left hand side of (|3"Tj) and 
the term on the right hand side in this expression. Moreover, we disregard the terms proportional to A and R" / R 
which become negligible as the flat direction starts evolving. In this limit, the masses of these two fields are 



m 2 



Higgs ^ j£2 



2 3£' 2 



+ 0(e±) , m L g =— +0( e 4) (32) 



where the terms in square parenthesis are of order e 2 . The time variation of these masses is adiabatic, and, therefore, 
does not lead to particle production. Also the mixing between the two modes is suppressed by higher powers of e , 
and therefore does not lead to any relevant physical effect. 

Finally, the spectrum of the fluctuations of this model also contains the two transverse polarizations of the gauge 
field, which are decoupled from the other modes (and from each other) at the linearized level. Their masses coincide 
with the longitudinal component given in (|32"f . 

Therefore, we conclude that a single flat direction does not experience nonperturbative decay, in agreement with 
what was argued in [111 ]. 



B. Multiple flat directions 

We now turn to the system with potential (|13p . It admits two flat directions, characterized by the background 
values 

(<M = (fo) = ^ > (&} = <<M = ^ 

As in the single flat direction case, we work in the phenomenologically relevant limit of 



e iS , (h) = <<M = iT5 eiS > <^> = (33) 



mR,mR,E',t', - F =, - ° = I < {F, G} (34) 
VF 2 + G 2 VF 2 + G 2 ) 

and we denote by e the ratio between any of the terms on the left hand side of (fM)) and any of the terms on the right 
hand side in this expression. 

The fluctuations 6<j>i of the scalar fields encode 8 real degrees of freedom. One of them is set to zero in the unitary 
gauge; the other seven degrees of freedom can be written in the form 

1 e lS 
54>i + 6(f> 2 = — (S r + iSi) , S(f>i - 5<f> 2 = — (S H + i F a) 
K H 



t>4 = ( '\ + ' <\ 



y(fe-*G«) (35) 
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We note that the same field a enters in 8<j>\ — 84>2 and 84>3 — (504. The difference between these two entries is precisely 
the Goldstone boson, which has been set to zero in the unitary gauge. 

The modes S r and 8i are the real and imaginary excitations of the first direction. They are decoupled from the 
other modes at the linearized level, and their mass matrix is identical to (J2U). Analogously, the modes 8 r and 8i are 
the real and imaginary excitations of the second direction. Their action is formally identical to that of 8 r and Si, 
upon the substitution of the background quantities and the model parameters related to the first direction with the 
corresponding quantities of the second direction. 

As in the previous case, we decompose the vector field into transverse and longitudinal components. The two 
transverse components are decoupled at the linearized level, and their mass is 



(F 2 + G 2 
4R 2 



(36) 



In summary, the above modes behave identically to those corresponding to the single flat direction; therefore, none 
of them are produced nonperturbatively. The action for the remaining modes, however, has no counterpart in the 
single flat direction case, and, as we now show, leads to nonperturbative production. Once the non-dynamical mode 
Aq is integrated out, we are left with four physical modes: the longitudinal vector polarization L, and the three 
fluctuations encoded in the second line of (|3"5|) . 

The two fields a and L are not canonical. In Appendix IBl we give the two linear combinations L 1; L 2 which are 
canonical in terms of a and L so that in matrix form, the action for the system is 



S^upicd = \ J d v d 3 k (A' 1 " A' + A'tjfA - A* if A' - Atfi 2 A) 



(37) 



where we defined A = (j>H > 8jj , L\ , , and where the matrices fl 2 and K are real and, respectively, symmetric 

and anti-symmetric. The exact expressions for these matrices are rather involved. In Appendix [B] we present them as 
an expansion series in e (defined in eq. Q34))). 

To compute the nonperturbative production of these modes, we need to find the eigenmasses of the physical modes, 
as well as the matrix T , defined in cq. (JT5J). For the eigenmasses, we find 



(F 2 + G 2 ) 

(F 2 m 2 + G 2 m 2 
F 2 + G 2 



(F 2 m 2 + G 2 m 2 ) 3 (W + G 2 ±> 



F 2 + G 2 R 2 (F 2 + G 2 f 



3(FG' - G F' 
R 2 (F 2 +G 2 f 



2 ZF 2 G 2 \iZ'-t' 



R? (F 2 + G 2 Y 



O e 4 



F 2 + G 2 
e 2 (F 2 + G 2 ) 



O e 4 



4i? 2 



0(e 4 ) 



(38) 



where terms outside of the parenthesis, and inside square parenthesis are, respectively, of zeroth and second order in 
e . We see that the system has two heavy and two light eigenstates. The fourth eigenmass coincides with that of the 
transverse vector modes (up to the accuracy of the present computation) and therefore we refer to the corresponding 
eigenstate as the physical longitudinal vector mode. The first eigenstate is instead the Higgs field of the model. As in 
the single flat direction case, these two fields have identical mass at leading order. 

For a single flat direction, the Higgs and the longitudinal vector mode were coupled only amongst each other in the 
action (|A15|) . Now, they are coupled to the additional light modes, which are absent in the single flat direction case. 
It is easy to verify that the eigenfrequencies of the light modes vary nonadiabatically with time. Namely, starting from 
the four eigenmasses (|38|) . we compute the comoving frequencies u>i — \J R 2 m 2 + k 2 , and we compute the adiabatic 
conditions lo[/ui 2 for momenta k comparable with the flat direction masses (hence, of order e in our notation). For 
the two heavy eigenstates, we find 
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This agrees with the general expression ([3]) valid for the heavy fields. However, for the light modes we find 



aJ 



0(1) 



(40) 



Hence, we see that the time variation of these eigenfrequencies is not suppressed in the limit of m -C ! As the 
numerical results presented in the next section show, this leads to the strong nonperturbative production of these 
modes. 

A second source of nonadiabaticity comes from the mixing in field space encoded in the (anti-symmetric) matrix T. 
We find, 



r 13 = - 



FG' — G F'f + F 2 G 2 (t,' - £') 
F 2 + G 2 

2 ) F 2 G 2 (S' - £' 



2 ^ 



>+0{e 2 ) 



F 2 S' + G 2 £' 
F 2 + G 2 



(FG' - G F'f + F 2 G 2 (£' - £' 



+ 0(e 2 ) 



Ti2 ; Ti4, r 2 4, r 34 — (e 2 ) 



(41) 



The terms within curly brackets are first order in e . As a consequence, two of the off-diagonal nonadiabaticity 
coefficients, defined in eq. ([26]) . are large 



^13 7 ^23 



0(1) 



(42) 



while the remaining ones are suppressed by an additional factor of e. This also leads to strong particle production (in 
particular, due to the A13 term, also the heavy Higgs mode is produced). 

We conclude this Section by noting that the physical longitudinal vector mode (more appropriately, the fourth 
eigenstate) has a negligible coupling with the other three modes of the system. This mode does not experience 
nonadiabatic production. 



V. NUMERICAL RESULTS 



We perform the numerical computation of particle production in the two flat direction case, with potential (|13[) . 

We focus on the system of fields |<5jr, &h, Li, L>2^ (as defined in the previous Section) for which nonperturbative 

production can be relevant. The background evolves as described in eqs. (|12p . with two additional equations for the 
amplitude and phase of the flat directions (identical in form as the first two of (fT2")l ). and with the potential of both 
flat directions appearing in the equation for the scale factor. Particle production is obtained from eqs. (| 18|) . 

The timescale for the evolution of the flat directions is set by their mass, so that terms of order e must be retained in 
the numerical equations solved. However, terms of higher order can be neglected. The coupling between the different 
modes is encoded in the matrix T, whose entries for the system we are studying are given in (|41|) . We then see that, 
if we neglect terms higher order in e, the longitudinal vector mode decouples from the other three. The mass of this 
mode varies adiabatically with time. Therefore, we can disregard this decoupled field in the computation, and we 
focus on the system of the three remaining modes. 

We should note that the potential (fl~3"|) is not bounded from below, due to the presence of the quartic terms 

proportional to A and A. Such terms dominate over the quadratic ones for large field values. In a complete model, 

we expect that higher order terms will be also present, and stabilize the potential. In our computation, we simply 

choose the parameters such that the quartic terms are subdominant (their presence is however crucial to generate the 

rotation of the flat directions in their complex planes). As we already discussed, the quartic terms can arise from 

.2 2 

one loop contributions to the effective potential and are assumed to have magnitude A ~ (and, analogously, for 

the second direction) in analogy to the baryon number violating operators in the Affleck-Dine mechanism where 
h 2 ~ 0.1 is some coupling constant. For definitcness, we set 
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This ensures that the quartic terms are initially subdominant, and that the two flat directions evolve towards the 
origin. 

We introduce the following dimcnsionlcss quantities (each denoted with an asterisk) , which are used in the numerical 
evolution: 



I 



e\$o\v 
m 



e|$o| 
F 



e|*o| 
G 



k 



Pip* 



e|$o| 

P4> 



A# = 



e 2 \$ \ 2 M? 



A* = 



72 



(44) 



The background equations, and the equations for particle production in dimcnsionlcss quantities are given, respectively, 
in eqs. (|D1|) and (|D2|) . One can see that the value of e factors out in these equations, and we only need to specify the 
value of |$o|/-Mp m the equation for the scale factor. 

We are free to set the initial value for the scale factor Rq = 1 , which implies 



*o 



1M 

l*ol 



(45) 



As we discussed in the introduction, flat directions can influence the thermal history of the universe only if they have 
a large initial amplitude. Therefore, we choose |$o|/A/ p = 10~ 2 in our numerical analysis. We also start with the flat 
directions at rest, an initial condition ascribed to Hubble friction. 8 At early times, the energy density of the universe 
is dominated by the inflaton contribution p^». As such, we start from the initial value p^* = 1 (see for example, the 
last expression of (|D1|) V As discussed in [llj, for a gravitational inflaton decay, the inflaton has yet to decay when the 
flat directions start evolving. At this stage, the inflaton is oscillating around the minimum of its potential, so that its 
energy density decreases as that of matter; therefore, we set p^* = 1/i? 3 in our numerical computations. 



(a) 



0.25 I — 
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- 
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-0.4 



(b) 
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Re 
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FIG. 1: The background evolution of the flat directions for the choice of masses m* = 10" , rh/m — 3.72 and initial phases 
So = 0.25 , So = 0.156. The initial values of the vevs are chosen such that iol/l'&o = 15 and | <B»o I = 10 _2 Af p . We show in 
the real and imaginary parts of the flat directions for the first three rotations of |$o| (solid red) and |$| (solid green); the right 
panel shows the same evolution as in the left panel, zoomed in on smaller values of the fields. 



For illustrative purposes, we show in figure[T]the evolution of the flat directions for some specific choice of parameters. 
We fix the soft masses to to* = 10 -6 , rh/m = 3.72 and numerically evolve the background equations (|D1[) with initial 
values .F*o = 2, G*o = 30 and So = 0.25, Sq = 0.156. In the left panel we show the evolution of the real and imaginary 



As clear from the parametrization J33}, this implies (F/R)' = S' = (G/R)' = £' = . The initial value for R' is obtained from the 
Friedmann equation. 
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parts of the flat direction vevs, normalized to their own initial values. In the right panel, we zoom in on the values 
closer to the origin, to better visualize the spiral motion of the two fields. 

We set the initial conditions for the Bogoliubov coefficients as cto = 1, /3o = 0, corresponding to an initial absence 
of particles. At the start of the background evolution, the masses of the light degrees of freedom are dominated by the 
tachyonic contribution i?"/i? 3 . This contribution rapidly becomes subdominant, and can be neglected when the flat 
directions start evolving. Our formalism for the production is valid only as long as the eigenfrequencies are positive. 
For this reason, we start the evolution of (fT8|) only once the eigenmasses have all become positive (we stress that the 
flat directions are still frozen at this stage). This amounts to neglecting the small gravitational particle production 
during the initial stages. 

From the numerical evolution, we obtain the occupation number densities rij of produced quanta, eq. (I22p . We can 
then compute the ratio between the energy density in these quanta, and that in the flat directions; in terms of the 
dimcnsionless quantities defined in (|4"4"|) , the ratio is 



7"prod 



Pprod 
Pflat 



167re 2 ) ( J ^* k* u i*(k*) n i(k*) 



x <R 




F„ R 



G* R 



R 2 F? 



< + f ^cos(4S) 



R 2 Gl 



^ 2 + T§ cos(4f]) 



(46) 



where the dimensionless eigenfrequencies are defined as lo^ = cji/ (e|$o|) ■ If '"prod becomes equal to one, we say 
that the flat directions have decayed. The present analysis is actually invalid at this point, since we are ignoring the 
backreaction of the produced quanta on the evolution of the flat directions (in particular, we do not account for the 
decrease in the amplitudes of the flat directions due to particle production). However, setting r pro d = 1 allows us to 
understand for which choice of parameters particle production is significant, and to estimate the decay time of the 
flat directions. It is useful to use as a measure of time the number of rotations of one of the two flat directions; for 
definiteness, we choose: N = (E — Eo) /2tt . We denote by iVdecay the value of N at which the ratio r pro( j equals to 
one. 

As we have mentioned, there are two different scales in this problem, set by the amplitude and the mass of the flat 
directions (the former controls the heavier eigenmass m\ in (|38p . while the latter the two light masses 7712,3). The 
timescale for particle production is set by the evolution of the flat directions, which in turn is governed by their mass; 
however, the intermediate matrices a and (3 evolve on a much quicker timescale, set by the largest eigenmass mi . 
Both scales need to be under control in the numerical simulations. The ratio between these two scales is the quantity 
e defined after eqs. (|34p. In realistic cases, these two scales differ by many orders of magnitude (GUT or Planck, 
vs. TeV scale). Simulations with such small values of e cannot be performed. Fortunately, as we discuss in appendix 
[El the occupation numbers of the two fields exhibit an approximate but accurate scaling with the masses and the 
amplitudes of the two directions. If we rescale both masses and both amplitudes by a constant factor 



{$o,$o} 
{m, m} 



{7$o,7*o}: 
{/im, nin} 



(47) 



then the occupation numbers scale as 



{ni , n 2 , n 3 } 



- rii , n 2 , n 3 
7 



Since e — > (p/7) e under (|47p . this is equivalent to saying that 



ni cx e 



772.3 oc e 







(48) 



(49) 



The numerical solutions obtained in the range 10~ 7 e < 10 -5 confirm this scaling behavior, as we show for a specific 
set of parameters in figure [2] We show the spectrum of the first (heavy) eigenstate for three distinct cases, which differ 
from each other only by the ratio between the masses and the amplitudes of the two directions (hence, only by the value 
of e). The three spectra coincide once n\ is rescaled according to eq. (|4*^]). and once the momentum is rescaled as k/e. 
Although we do not show it here, we also verified that the occupation numbers for the light eigenstates do not change 
with e. The system of equations for particle production is too involved to prove this scaling analytically. However, 
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FIG. 2: Occupation number for the first (heavy) eigenstate after three rotations of the first direction, for a specific set of 
parameters: |<E> | = 10 _3 M P , |$ |/|$o| = 4, rh/m = 3.72, E = 0.25, E = 0.156. The larger figure shows the occupation 
number as a function of the momentum (rescaled as in eq. (|44p ) for three specific choices of the (rescaled) flat direction mass 
m„. The insert shows the same result, but plotting ni/ro* vs. fe«/m». In terms of these variables, the three curves overlap. 
Since m» oc e, this confirms the scaling (|49[l for the occupation number, and the fact that the momenta of the quanta produced 
also scale as e . 



in appendix [El we see that the scaling "emerges" from the equations under a few assumptions (|eay| <C \u>i oc\j\ and 
|e/3y| <§C \u>i for j = 1, 2, 3) which are verified numerically. 

Perhaps the most convincing argument in favor of the scaling (|49|) is the fact that the nonadiabaticity parameters for 
the light states, aj 2 /wf and uj' 3 /u}^, are not suppressed in the limit of small e , cf. eqs. (|4"0)) . The analogous parameter 
for the heavy state is suppressed. However, the production of this state takes place through the coupling with the two 
lighter ones (cf. eq. (|4"2")l . Since u>i oc e , u>2,3 oc e, equipartition of energy then suggests that n\ is suppressed by e 
with respect to 72.2,3 ■ 

Due to these scaling properties, the ratio between the energy density of the produced quanta and that of the 
background flat directions scales as (see appendix IE]) 



r-p r od (50) 

T 



This implies that r pro( j is of the form 



m m 




Vod * T^r-TZ- X / -, N] (51) 



The multiplying function / can be computed numerically for cases in which the amplitudes and the masses of the flat 
directions differ only of a few orders of magnitude (the approximate scaling holds only as long as e < 1). The scaling 
(j5"Tj) then gives the production for realistic cases (e as small as 10~ 15 ) for which numerical evolution is not feasible. 

The function / obtained numerically exhibits a strong dependence on the parameters. The rate of particle production 
grows exponentially within a range of parameters, while it is negligible otherwise. Therefore, eq. (|5ip can be cast in 
the form 

r prod ^C-J^-W^, (52) 
|$o| |$o| 

where C and a are two time-independent quantities that are functions only of the ratios rh/m and |$o|/|^o|- The 
growth rate, a, rapidly approaches zero outside the range for which particle production takes place. 
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In figure [3J we show the numerical values for a obtained from the numerical simulations for three given mass ratios, 
and for a range of ratios between the two amplitudes. 




FIG. 3: The dependence of the growth exponent a defined in (|52|l on the initial ratio of the vevs |<&o|/]$o], for three different 
mass ratios. 



As is clear from (|52j) . iVdocay is inversely proportional to a . 



^decay - - log 10 -fTZ • (53) 

tr \ C mm J 

We show this quantity in figure 31 for the same ratios of masses and amplitudes used in figure [3J The overall mass 
scale is set to m/e = 10 4 GeV (the mass of the first direction). The overall scale of the amplitudes is fixed by rescaling 

/ ~ \ 1/2 

the results to (|$ | |$o|J = l(T 2 M p . 

We see that the production of particles is extremely fast, provided the initial amplitudes of the two flat directions 
are not too different. Not surprisingly, a flat direction with a significantly smaller amplitude can be neglected, so that 
one is effectively back to the single flat direction case, for which we know that particle production is absent. This 
can be seen more clearly in figure O where we show the sum of the number densities of the produced quanta (for a 
given value of the momentum) in the upper panel, and the evolution of the amplitudes of the flat directions in the 
lower panel. We note that particle production takes place in a stepwise manner, whenever the amplitudes of the two 
directions are comparable to each other. The overall growth is exponential with time, as parametrized in equation 
([12]), once we average over complete rotations (this is also how the quantity <r, plotted in figure [3] has been obtained). 

The fact that particle production takes place whenever the two directions have comparable amplitude does not mean 
that the two amplitudes need to be equal (or nearly equal) initially. This is due to the fact that the orbits of the two 
flat directions are strongly elliptical (cf. figure [1]), so that each amplitude varies significantly during its rotation. As a 
result, we see in figure 2] that production takes place for a significantly wide range (about four orders of magnitudes) 
of the ratio of the initial amplitudes. We actually see that the central point of this interval moves to larger values 
of |$o|/|^o| as m/m increases. To understand this, imagine that the second direction is heavier than the first one. 
Then, there is an interval of time during which the first direction is frozen, while the second one evolves (these are 
the times for which m < H (t) < m). During this interval, the amplitude of the second direction (more appropriately, 
the maxima and the minima of the amplitude) decreases, due to Hubble friction, as i?~ 3 / 2 , while the amplitude of 
the first direction remains constant. Therefore, to have maximal overlapping between the two amplitudes when both 
directions are evolving, the second direction should start with a higher amplitude than the first one. The value of 
l^o |/|$o | leading to maximal overlapping (and maximal production) increases as m/m increases, as figure H] clearly 
shows. 
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FIG. 4: Number of rotations of the first vev where the production criterion r pro d = 1 is satisfied and its dependence on the 
initial ratio of the vevs |<&o|/|^o|) for three different mass ratios. The numerical analysis was made for the first twenty rotations 

of the first fiat direction and the results were rescaled to m/e = 10 4 GeV and J |<I>oj |$o| = 10~ 2 M P . 



Finally, we observe that the large production of particles can be traced to the adiabatic parameters. The leading 
order adiabatic matrix elements A22, ^33, A13 and A23 are all of order one in the series approximation (|34p and 
generically are greater than one during much of the evolution. We show the root mean square (RMS) of the leading 
order adiabatic matrix elements in figure [6l for the mass ratio rh = 3.72 m. As explained in Sections HT1 and IIII1 a 
necessary (but not sufficient) condition for the production of quanta is realized when any element of the adiabatic 
matrix is greater than one. By comparing figure U with figure [SI we see that the results support this assertion. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper, we have computed the nonperturbative decay of systems of multiple flat directions. Our results agree 
with the conclusions found in [ll[. However, we improve over that study by including gauge fields in the computation. 
We worked out in detail the production of particles for systems with one or two U(l) flat directions. Specifically, we 
studied the linearized theory of fluctuations around such backgrounds, and we computed the production of particles 
due to the time dependency of the physical eigenstates and eigenmasses of the systems. We also showed that this 
computation can be readily applied to SU(N) directions, as long as the gauge fields do not develop background 
expectation values. 

Technical difficulties associated with this study revolved around the fact that one needs to find the spectrum of the 
fluctuations around a time dependent background. Firstly, one finds that, even in the unitary gauge, the longitudinal 
component of the gauge boson remains coupled to the "matter" fluctuations (the Higgs and light fields) at the linearized 
level. This coupling was not included in the analogous computations of [2lT[2^|. which adopt a gauge choice equivalent 
to that chosen here. However, this coupling needed to be included for consistency, and forced us to perform some field 
redefinitions (see Appendix [Bl for details). 

Secondly, the quadratic action for the perturbations contains nonstandard quadratic kinetic terms, and "mixed 
kinetic terms" of the type K^S^ Stfi'j + X^Sfy S(j)j , where the matrixes K and X depend on the background and are 
therefore time dependent. We provide a general formalism to diagonalize (and eventually quantize) such system, which 
is a necessary step to compute particle production. In particular, we show how particle production can be obtained 
without explicitly performing the rotation that eliminates the mixed term proportional to X. 

Thirdly, the resulting system of equations for particle production is too involved to be solved analytically. In 
the case of two flat directions, the relevant computation can be reduced to a complicated system of four coupled 
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FIG. 5: Upper panel: Number density of the produced quanta for a given momentum fc* = 10~ . Lower panel: amplitudes of 
the two fiat direction. We note that particle production occurs whenever the two amplitudes are comparable, as discussed in 
the main text. The parameters chosen for this evolution are m* = 10 -6 , F, — 2, fh/m = 3.72, |3 > o|/| ( I > o| = 15 



fields. To compute particle production, one needs to know the eigenvectors and the eigenmasses of this system. The 
analytic computation of these quantities is unpractical. On the other hand, the system possesses two mass scales, 
one given by the masses m of the flat directions, and one by their amplitudes |$| . While m is expected to be of 
the order of the electroweak scale, |<E>| can be as large as the GUT or the Planck scale (indeed, very large vevs are 
necessary for the flat directions to play an important role in the thermal history of the universe jll|). Any numerical 
evaluation of particle production needs to control both scales; this is also an untreatable problem, for realistic values 
of e = m/|$| . Fortunately, the computation of the eigenmasses and eigenvectors can be performed analytically as an 
expansion series in e . We could then perform a numerical evaluation of the resulting formulae for intermediate values 
1CU 7 < e < 1CU 5 . The results of these evaluations showed that the particle production exhibits a clear scaling with 
e (we also presented several analytic arguments in support of the observed scaling). The most important outcome is 
that the non adiabaticity in the evolution of the eigenmasses and the eigenvectors is not suppressed in the limit of 
e — > , as was the case for the toy models considered in [l(| HH, [H| ■ As a consequence, the flat directions can undergo 
very fast nonperturbative decay even in this limit. 

The results of [TTf| for the non-perturbative decay of flat directions have been discussed in some recent papers 
[20l. l2ll I22I I23I . Ref. [2(| agrees with [ll[ that the effect is absent when only a single flat direction is excited; however, 
contrary to it was argued there that the nonperturbative decay is unlikely in the case of multiple flat directions. 
The main argument against this is that, according to [20], (i) flat directions present hierarchical vevs, and (ii) there 
is no nonperturbative particle production whenever the vevs of the two directions differ by more than one order of 
magnitude (for more details, see their eq. (27), and the discussion immediately afterwards). The claim (i) is based on 
the fact that the higher order terms 4> n /M n ~ 3 that we have discussed in the Introduction arise at different values of 
n for different directions. As we already mentioned, the fact that one such term is allowed does not necessarily mean 
that it is present. On the contrary, we have already seen that, for the simplest inflaton model we are considering, 
all terms with n < 11 must be forbidden. Therefore, either such terms are present (in which case the vevs of the 
flat directions are too small), or are forbidden (in which case we cannot use this argument to see at which order any 
given direction is lifted). Also the claim (ii) is not substantiated by numerical computations, or by a clear proof. Our 
explicit computations show that this conclusion is actually incorrect, since we find particle production for a range of 
four orders of magnitude among the initial values of the two amplitudes. 

In contrast, ref. 



and agrees with [11 



2l| also finds that the effect of particle production can be present in the multiple flat directions case 
. However, it is argued there that particle production can be present also for a single flat direction. 



This conclusion is based on a specific example (their Section IV), where nonperturbative production takes place. Ref. 
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FIG. 6: The root mean square (over physical time) of the four leading adiabatic matrix elements A22, A33, A13 and A23 of the 
modes with momentum k = 0, for five rotations of the flat direction $. While the RMS values are not directly correlated with 
the growth exponent a, they indicate the parameter space where particle production may be possible. 



[2l| claims that there is one flat direction in this example. We agree that there is nonperturbative production in this 
model. However, this model contains two flat directions. This is clear from the fact that the background fields in 
this model have two physical (i.e., that cannot be removed by a gauge transformation) phases, their a and 7, that 
evolve independently with time (we count the number of flat directions, based on the number of independent physical 
background fields involved). 

Ref. [22| studied whether realistic MSSM flat directions can undergo nonperturbative decay. It was concluded there 
that this effect does not take place for a single flat direction LLe c or u c d c d c . It was also shown that nonperturbative 
decay can take place for the simultaneous presence of the directions LLe c and QLd c , but not for the simultaneous 
presence of LLe c and u°d°d c . These conclusions agree with those of [ll[. In particular, it is easy to understand why 
the nonperturbative decay does not take place when the LLe c and u c d c d c directions are excited. In the presence of 
these two directions, the D— term potential involves 28 fields (the excitations of the fields which acquire nonvanishing 
background values, plus those related tho these fields by off-diagonal gauge generators). These two directions break all 
the SM gauge symmetries. As a consequence, there are 12 goldstone bosons, and 12 Higgses in the spectrum. These 24 
fields, plus the 4 fields representing the real and imaginary excitations of the two flat directions, exhaust all the fields 
in the spectrum. Therefore, there are no additional light fields in which the heavy fields can rotate into. This is the 
reason why these two directions alone were not considered in [ll[ ; however, these two directions allow for the presence 
of a third one, of the QLd c type. In this case, there are sufficient fields in the spectrum to make nonperturbative 
decay possible. 

Finally, ref. [23j | claims that flat directions cannot undergo nonperturbative decay, due to some charge conservation. 
Ultimately, the conserved charge is the angular momentum of the background flat directions (in field space; for many 
directions this is related to the baryon or lepton number). We do not claim that the D— term interactions remove 
this charge. However, they can redistribute it from the homogeneous condensate (k — 0) to the quanta of the fields 
coupled to the condensate through the D— term potential. This effect is missed if one only considers the homogenous 
fields. 

To conclude, we confirmed the results of that two or more flat directions possess one nonperturbative instability 
that can result in a much faster decay than what one would simply argue from their perturbative interactions. As it is 
typical for preheating studies, this instability is due to the nonadiabatic evolution of the eigenmasses and eigenvectors 
of the spectrum of the theory around the background flat directions. Once the production becomes significant, the 
nonlinear interactions of the produced quanta, and their backreaction on the homogeneous flat directions (besides the 
obvious decrease of the energy of the flat directions) will likely become important. This is a nontrivial issue, due to 
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the complication of the MSSM potential [11], which deserves further study. 
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APPENDIX A: ONE U(l) FLAT DIRECTION 

In this appendix we study the flat direction composed of two scalar fields <pi , 4>2 with the Lagrangian ([T]) and the 
potential Let us first comment on the simplest case of a U(l) symmetric theory with a single complex field x 
(giving rise to the standard Higgs mechanism). Denoting, respectively, by v and £ the amplitude and phase of this 
field, and taking the covariant derivative as in |T]), with q x = 1 , one finds the following coupling between the scalar 
and the gauge fields 

|2 „,2 ^ - - e 



\D X \ A DV[d li Z--A ll ) (Al) 

It is then conventional to choose the unitary gauge in which the phase £ is set to zero. In this gauge, the linear coupling 
between the vector and the scalar is absent. If we then expand v in a (nonvanishing) vacuum expectation value plus 
fluctuation, the spectrum of the theory consists of a Higgs field, and a massive gauge boson, which are decoupled from 
each other at the quadratic level. 

Let us now turn to the case of our interest. We decompose the two fields entering in the potential © as 

j. f M0+a) J. 9 



where /, g, a and (5 are real fields. The parametrization chosen for the phases is related to the fact that a is the only 
quantity changing under the U(l) transformation (since the two fields have opposite U(l) charge: q\ = 1, q 2 = — 1). 
We can completely fix the gauge freedom by going to the unitary gauge, in which a is set to zero. The coupling 
between the gauge and the scalar fields now reads 



\D<k | 2 + \D<f) 2 \ 2 D — 2 {-f + g 2 ) d»P A, + { f+g 2 ) A^ (A2) 



so that in general a linear term in A^ is present even in the unitary gauge. 
We decompose the scalar fields into background values plus fluctuations, 



(A3) 



and we restrict our attention to the flat direction background 

{f) = {g) = ^= , </?)^S (0 1 ) = (^) = Z-e iS (A4) 

(for which the D-term potential vanishes). Since the amplitudes of the two complex fields are equal, the linear term 
in Afj, in (|A2|1 vanishes on this background, so that we can consistently set = . 
From these background assignments, we get the background action 



Sbck =2 / 



F' 2 + F 2 S' 2 - ( m 2 R 2 -?pjF 2 - ^F A cos(4E) 



(A5) 



from which the first two equations in (|12p follow. 

We next study the linearized theory for the fluctuations. This amounts to expanding the total action to quadratic 
order in the fluctuations 5<f>i and in A^ (since the gauge field has vanishing expectation value, it is treated as a 
fluctuation). The perturbations of the scalar fields contain three real modes in the unitary gauge. It is convenient to 
redefine them as 

'^t ■ s ' JsT ~7T ■ »-j <-> 
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in terms of which 



e iS 1 i 1 

ii + <5</>2 = (r + i a) = — (cos Er - sin £ <r) H (sin Er + cos £ cr) = — (<5 r + i Si) 

R R R R 

i>l - ^02 = <5tf -fr- 
it 



(A7) 



It is also convenient to decompose the spatial components of the vector field in a transverse plus longitudinal part 

Ai = Af + d % L (A8) 

where diAj = . 

As we now show, these combinations have an immediate interpretation. Indeed, the quadratic action for the 
fluctuations splits in three separate pieces 



S {2) = S flat [8 r , Si] + S X [Af] + S|| Higgs [6 H , L, A q ] 



(A9) 



The modes <5 r , Si can be interpreted as the real and imaginary excitations of the flat direction. Their quadratic 
action is 



Sflat = - / d 4 x 



(<9,A) 2 + (d^i) 2 - (S r , Si) 



m 2 R 2 - ^ + ^ F 2 cos (2£) 
-^F 2 sin(2S) 



-^F 2 sin(2£) 



(A10) 



The quadratic action for the transverse vector modes gets its contributions from the kinetic gauge term —F 2 /4 and 
from the second term in (IA2D. We find 



Si 



1 



d A x 



A-i A i 



e 2 F 2 



-Ai A i 



(AH) 



where prime denotes (conformal) time derivative. 

The remaining fluctuations enter in the quadratic action 



I 



IIHiggs 



d 4 x S' H S' H - (d t 8 H ) (diS H ) + (diL') (diL') + {diAo) (^Ao) - 2 (d^o) (d t L r ) 

e 2 F 2 e 2 F 2 

-2 e F £' A S H -— (diL) (d t L) + ——A n A 



e 2 F 2 Q „., II" . .,.> A 



m 2 R 2 - — - S' 2 - —F 2 cos(4E) ) S H 5 h 



(A12) 



The second line is obtained from the coupling term (|A2[) . We see that the only linear term in the gauge field couples 
Aq to the mode Sr ■ We also see that the mode Aq is nondynamical, and can therefore be integrated out. After Fourier 
transforming the spatial coordinates, the equation of motion for Aq gives 



An =4 



eFT,' S H + k 2 L 
4 k 2 + e 2 F 2 



(A13) 



We can substitute this solution back into the Fourier transform of (|A12[) . In this way we obtain the action for the 
two dynamical modes Sh and L . The field L is not canonically normalized. The canonical variable is 



-4 ii 



ekF 



Ve 2 F 2 + 4 k 2 



L 



In matrix notation, the action takes the form 

S'llHiggs = 



X - J d v d 3 k (A't A' + A'^KA - A^KA' - A t fi 2 A) , A = 



Sh 



(A14) 



(A15) 
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where the matrix K is antisymmetric, Q 2 is symmetric, and 
2kT,' 



K12 = 



(tt 2 ) 22 

(« 2 )l2 



Ve 2 F 2 + 4 k 2 
e 2 F 2 



4 

e 2 F 2 



f- TO i? — 

4 fc 2 



i?" AF 2 



e 2 F 2 + 4 fc 2 



i? 2 

e 2 F 2 



kF 2 



Ve 2 F 2 + 4 fc 2 



6e 2 £' 



cos(4S) - 
/c 2 + m 2 R 2 



3e 2 F 2 — 4 fc 2 
e 2 F 2 + 4 fc 2 
AF 2 



cos(4 E) - E' ; 



2 ci/2 



3e 2 F 



e 2 F 2 + 4 k 2 



EL 

~F 



A sin(4S) 



(A16) 



To eliminate the "mixed kinetic terms" (those proportional to K) we perform a field redefinition 'J' = 7?. A, where 
7?. is orthogonal (so that the quadratic kinetic term remains canonical) . Under such a transformation 



i^'t^r' = I ( A 't A' + A't R T R' A - A*i? T R'A' + A^ R' T R' A) 



(A17) 



where in the third term we have used the fact that R T R 1 is antisymmetric (due to the orthogonality of R) . If we now 
identify 

TZ T W = K 

which gives K T K = 7Z' T 7Z' , the action (|A12|) becomes 



<S||Hig g s = \j dr l 63 k [*' f *' _ * f 11 ( Q2 + kTr ) nT *] 



To proceed, we need to diagonalize the frequency matrix 1Z (il 2 + K T K^j 1Z T : 

c T n{n 2 + K T K)n T c = u J 2 ._, - r = c T c< 



(Ai8) 

(A19) 
(A20) 



diagonal 

As outlined in Section HVl both the eigenfrequencies and the matrix T are necessary for computing the production of 
quanta from the evolving flat direction. The matrix TZ is specified by the condition (|A18[) and by the orthogonality 
requirement. However, the explicit knowledge of TZ is not necessary to solve the eigenvalue problem (IA20[) . Firstly, 
since TZ is orthogonal, the eigenvalues of 1Z (ft 2 + K T K^j 1Z T coincide with those of f2 2 + K T K. Secondly, if we simply 
diagonalize the matrix SI 2 + K T K, 

£ T (fl 2 + K T K) £ = diagonal 



(A21) 



we can write 



C = 1Z£, => T = C T C =£ T +£ T K£ (A22) 
where the condition (|A18|) has been used. Therefore, we simply need to know the eigenvalues and eigenvectors of 

n 2 + k t k. 

The problem is actually simplified by noting that the off diagonal terms are suppressed in the phenomenologically 
relevant cases, namely when the conditions pip are met. Disregarding the terms proportional to A and R" (which 
become negligible as the flat directions start evolving), we find 



s> 2 



K T K = 



[k 2 



m 2 i? 2 + 3£' 2 ] 



0( e 4 ) 



0( £ 3 ) 
^l+\k 2 ]+0^) 



(A23) 



where the terms in square parenthesis are of order e 2 . The diagonal elements (the eigenfrequencies of the system) 
vary adiabatically, u>' /oj 2 = O (e) . 

We can also expand the matrix T obtained in (|A22|) . 



r 12 
-r 12 



r. 



Akm 2 R 2 £' 
eF (m 2 i? 2 + 3£' 2 ) 



3F' 



F (m 2 i? 2 + 3E' 2 ) 





R'\\ 


(4- 





0(e 4 ) 



(A24) 



so that we see that the dominant term is of order e 2 . As a consequence, the matrices /, J entering in the equations 
(fT5|) for the particle production are 

0(e 2 ) 
0(e 2 ) 



I = 



J = 



C(e 4 ) 
C(e 4 ) 



(A25) 
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APPENDIX B: TWO U(l) FLAT DIRECTIONS 



Now consider the two flat directions composed of four scalar fields with the potential (jT3|> . This potential can 
accommodate two flat directions, characterized by the background assignments 



2R 



G 







(Bl) 



We parametrize the fluctuations as follows 



F + r + 5 H 



i (S+cr/F+a) 



2R 



2R 



lH_ e i(-£+a/G+a) 



F + r-S H Us+a/F-g) 

2R 

G + f-8 H j (S+a/G-5) 

2R 



(B2) 



where all the fields are real. The scalar fields cj>i and <p^ have positive (+1) U(l) charge, while (f>2 and </>4 have negative 
(— f ) charge. For this reason, the only linear combination of the phases in (|B2I) which is gauge variant is a + a which 
we set to zero in the unitary gauge. We also define a — (a — a)/2. The background value of the latter combination is 
forced to zero by the equations of motion. 

Calculating the background action yields two copies of (|A5jl . As in the single flat direction case, we decompose the 
transverse and longitudinal parts of the vector field as: 



A, = Af + d t L , 



(B3) 



where diAj — 0. At the linearized level, the transverse part of the vector field forms a decoupled system with action 



AJ' AJ' - ( di Aj) ( di Aj) 



e 2 (F 2 + G 2 ) rp rp 

Z *H A i 



As in the previous case, the two combinations 



e f 

h + &4>2 = (r + i a) —— = — (S r + i Si) 



(B4) 



(B5) 



are the real and imaginary excitations of the first flat direction. The quadratic action for these two modes is identical 
to (|A10p . Analogously, the two combinations 



94 = 



e iS f /- ~ \ 
(f + ia) — = — [S r + iSij 



(B6) 



are the real and imaginary excitations of the second direction. Their quadratic actions is formally of the type (j A10[) 
with all the quantities referring to the first direction substituted by the analogous ones of the second direction. 
The remaining 5 modes are coupled into the quadratic action 



S, 



coupled 



P 2 (p2 , n2\ 

+ (F 2 + G 2 )(d^a) (d»a) 1 ^ - ' (diL) (diL) 

-2 (diA ) (diL') - 2 e (f S' S h + G t' S H ) A 

-e (F 2 - G 2 ) [A a' - (diL) (did)] + 4 (f E' S h - G £' 6 H ) «' 

e 2 (F 2 + G 2 ) e 2 / - \ 2 

+ - 1 / > A A ~ e -(FS H + G S H ) 

R" A \ 

m 2 R 2 - — - E' 2 - -F 2 cos(4E) J 5 H 5 H 

- L 2 R 2 -^- E' 2 - ^G 2 cos(4E) N ) 8 H S H 



(B7) 
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To proceed, we Fourier transform this action with respect to the spatial coordinates, and we integrate out the 
nondynamical field Aq . We find 



A = 



1 



fe 2 + f (F 2 + G 2 ) 



(F z -G 2 ) a' + k 2 L' + e 



(B8) 



We then insert the solution (|B8[) back into the Fourier transformed action. The kinetic term of the resulting action is 



Skin = ^ I dr]d 3 k 



l^| 2 + l^| 2 + 



e 2 k 2 (F 2 + G 2 ) 2 f 4 [e 2 F 2 G 2 + k 2 (F 2 + G 2 )} ^ 



e 2 {F 2 + G 2 )+Ak 2 
2ek 2 (F 2 ~G 2 ) 



e 2 {F 2 + G 2 )+4k 2 



(L a + L a ) 



e 2 (F 2 + G 2 )+4k 2 

We canonically normalize the a and L modes through a series of redefinitions. We first define 



L = 



e k VF 2 + G 2 
yje 2 {F 2 + G 2 )+4k 2 

in terms of which the kinetic term reads 

f 



L 



a=2, 



Ski 



kin 



drjd^k 



l^| 2 + l^| 2 + |i'| 2 + |af 



e 2 F 2 G 2 + k 2 (F 2 + G 2 ) 
I e 2 (F 2 + G 2 ) +4fc 2 

k (F 2 - G 2 ) 



■ a , 



yJ{F 2 + G 2 ) [e 2 F 2 G 2 + k 2 (F 2 + G 2 )} 



L a + L a 



The mixed double derivative is eliminated by 



Li 



L — a 

~7T ' 



F 2 



L + a 

~7T 



Finally, the canonical variables are 



^(F 2 + G 2 ) (e 2 F 2 G 2 + k 2 (F 2 + G 2 )) + k (F 2 - G 2 



1/2 



[(F 2 + G 2 ) (e 2 F 2 G 2 + fc 2 (F 2 + G 2 ))} 



1/4 



L 2 = 



^{F 2 + G 2 ) (e 2 F 2 G 2 + k 2 (F 2 + G 2 )) - k (F 2 - G 2 



1/2 



[(F 2 + G 2 ) (e 2 F 2 G 2 + k 2 (F 2 + G 2 ))] 1/4 
In terms of these fields, the action is formally of the type 

S^pied = lj dr l d3k ( A ' f A ' + A ' tKA - A f FA' - AW A) , A 



fs H \ 

~5h 
Li 
\L 2 J 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 



where the matrices fi and K are real and, respectively, symmetric and anti-symmetric. 

The exact expressions for these matrices are rather involved. However, we can present them as an expansion series 
in e (defined in eq. (f3"4")) ; as in the single flat direction case, we can neglect the terms proportional to A, A and R"). 
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Wc find 
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K 23 
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f Gg 
GE' 



+ 



\/2\/F 2 + G 2 
Ft' 



fcE' 



fcE' 



3F 2 + G 2 
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3F 2 + G 2 



( V2VF 2 + G 2 

J \ 

\V2VF 2 + G 2 J 

fc (F 2 - G 2 ) 
2eFG (F 2 + G 2 ) 3/2 



2V2eF V ^ 2 + G 2 
fcE' /F 2 + 3G 2 



2V2eG V ^ 2 + G 2 

fcE' / F 2 + 3G 2 
F 2 + G 2 

G' 



2V2eG 



G2 ? + ^G 



+ 0(e 3 ), 
+ 0(e 3 ), 

+ 0(e 3 ), 

+ 0(e 3 ), 



(« 2 )ii 
(« 2 )i2 



(f2 2 )i3 



(0 2 ) 14 



e 2 F 2 



e 2 F G 



Jr + m 2 i? 2 + f 3 iT f 2 ) E' 2 



F 2 + G 2 



AFG 
F 2 + G 2 



E' E' 



G(e 4 ), 



3F 2 GE' 



V2 (F 2 +G 2 ) 3/2 K~F ~G 
3F 2 GE' / F' G' 



V2 (F 2 + G 2 ) 3/2 G 



+ 0(e 3 ) , 
+ 0(e 3 ). 



(n 2 ) 22 



(« 2 ) 23 



(0 2 ) 24 



(^ 2 ) 33 



e 2 G 2 



■ + 



fc 2 +m 2 fi 2 -[ f l~ 3 g jE' 2 



F 2 + G 2 



+ 0(^ 4 ), 



3FG 2 E' 



(F 2 + G 2 f 



F' G' 



3/2 \ F Q 



0(e 3 ) 



3FG 2 E' 



F' G' 



V2 (F 2 + G 2 ) 3/2 \ F GJ\ ° {e3) 
e 2 (F 2 + G 2 ) | j ek{F 2 -G 2 )VF 2 + G 2 \ 



8FG 



fc 2 + 



(« 2 ) 34 



+ 

e 2 (F 2 + G 2 ) 



F 2 (m 2 i? 2 - E' 2 ) + G 2 (m 2 i? 2 - E' 2 ) 3 (GF'-F G') 2 



2(F 2 + G 2 ) + 2(F 2 + G 2 ) 2 



+ G(e 3 ). 



fc 2 (F 2 -G 2 ) 2 | F 2 (m 2 i? 2 -E- 2 )+G 2 (m 2 i? 2 -E- 2 ) | 3(Gf /_ FG /) 
"""" + 2 (F 2 + G 2 ) 



16 F 2 G 2 



2 (F 2 + G 2 ) 2 



(« 2 ) 44 



: (F 2 + G 2 ) / efc(F 2 -G 2 )VF 2 + G 2 | 



8FG 



+ 



. F 2 (^^-E^+G 2 ( TO 2 i? 2 -E' 2 ) 3 (GF'-FGT 
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+ 0(e 3 ) 



(B15) 
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In these expressions, terms outside parenthesis, within curly parenthesis, and within square parenthesis are, respec- 
tively, of zeroth, first, and second order in e . 

We then proceed as in the previous Appendix by computing the eigenfrequcncies, and the matrix T. The leading 
expressions for these quantities can be found in the main text. 



APPENDIX C: SU(2) FLAT DIRECTIONS 

The linearized computation done in the previous two appendices for the U(l) case can be easily extended to 
nonabelian flat directions of physical interest. The reason for this is that the differences due to the nonabelian 
structure arise only at higher than quadratic order in the gauge fields. Such interactions are neglected in the linearized 
study of fluctuations, as long as the gauge fields involved do not have any background expectation value. Therefore, 
in the linearized study, the systems of physical fields around nonabelian directions can be treated as a series of copies 
(one per generator of the group) of those encountered in the abelian case. We show how this is realized in specific 
examples involving one or two SU(2) flat directions. 

Let us start from the case of a single flat direction. We introduce two SU(2) doublets, with the potential (O, and 
with the background values 

F je fi\ f JE (a 



2R^[0) ' <« = 2« e 'HlJ (C1> 

breaking all the SU(2) symmetries. We also assume that the three gauge fields have vanishing background expec- 
tation value. There are eight real fluctuations of these two complex doublets. It is convenient to parametrize them 
as 

_ F + r + S H jE / 1 + i \ _ F + r-S H iE / ^ + * (X - 



2R ° \ L&ii+i (s 1+ i~S 2 ) I ' 02 2i? C I ! + 



(C2) 



Under the infinitesimal SU(2) transformations & —> U (pi, with U = l + ia a a a /2 (where a a are the Pauli matrices), 
only the modes Si are gauge variant, transforming as 

Si -> 5i + ^ . (C3) 

In the following, we fix the gauge completely by setting Si = . 

The background action is the same as in the U(l) case given by eq (|A5[) . For the fluctuations, we redefine the fields 
r and £3 as we did for r and a in the first line of (|A7[) . We then see that the action can be split in several pieces, 
analogous to those found for the U(l) case, 

3 

S = S l]mess (Al L\ S x ) + S miggs (Al L 2 , S 2 ) + S miggs (A 3 , L 3 ,S H ) + 5 flat (5 r , 6 t ) + ^ Sj_(Af T ) , (C4) 

a=l 

where Sa^t is given in eq. (| A10|l . the S± actions are three copies of (|A11[) . and the SiiHiggs actions are three copies of 
(|A12p . This confirms what we argued in the main text, namely that the coupled system for the SU(2) flat direction 
consists of three copies (one per generator) of that obtained in the U(l) case. 

This is true also for multiple flat directions. Consider the case of four complex doublets with the potential ([15)) . 
Specifically, we are interested in the case in which two flat directions are present. The doublets </>i and <f>2 have 
background values as in (|C1|) . The doublets $3 and (f>4 have analogous values, with F replaced by G and S by S. 
Similarly to what was done in (|C2|) . it is convenient to parametrize the 16 real fluctuations of these fields as 



F + r + 5 H JS / l + i (#+4) 
2R C \ tiMh. +i (s 1 + il 



G + f + S H _ it ( l + i (^+j 3 ) } ^ G + ?-S H et ( + l (71 -17a) 




2R V 2i ^ + * (7i+*72) )' ^ 2R U + Hi-^j 



(C5) 
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In this case, the gauge dependent combinations of the perturbations are Si + ji. We fix the gauge completely by 
setting them to zero, and we define a, = (Si — 7i)/2. 

The background action yields two copies of (|A5p . Transforming the fields (r,S 3 ) and (7,73) as in the first line of 
(|A7[) . we can express the quadratic action in terms of the separated actions obtained in the U(l) case, 

S = 5'couplcd(^4o, L , Si, 71, Oi) + <S CO upled(-^0> ^ > ^2, 72, a 2) + ^coupled (Aq , L 3 , Sh, Sh, ^3) 

3 

+ Sflat(<5rA)+Sfiat(^A)+E S±(^ T ), (C6) 

a=l 

Here as well, the coupled system also consists of three copies of the coupled system obtained in U(l) case. 



APPENDIX D: EQUATIONS IN PROGRAM UNITS 



We list here the closed set of equations that we solve numerically in Section [V] In terms of the dimensionless 
quantities (J44J) , the background equations are 



R 



A* 



K + ( ™lR 2 - - S' 2 ) F* + ^Fi i cos(4S) = 



(i?E')'-^#sin(4E)=0 



G" 



fhlR 2 - 



^-£>A G. + ^G» cos(4E)=0 



( G 2^y_^l G 4 sin ( 4 £ 



•12 



R^_ R_ 

R + R 2 



■ 4-7T 



M p 2 



\ p4 
_2 i^2 1 A * 



cos (4E) + m 2 Gl + cos ( 4^ 



where prime now denotes derivatives with respect to rj* 
Eqs. (fl8|) for the produced quanta read 



a' = (—ioJ* — /*) a + f — J*j P 
f¥ = (tw„-/0)9+ {i;- J*) a 
where we have introduced the dimensionless frequencies 



= diag (wi*, w 2 *, w 3 *) , 



-fc 2 



i? 2 



(Dl) 



(D2) 



(D3) 



with the eigenmasses (|38|) , and the dimensionless matrices 



_ii 



e l$r 



(D4) 



with the r matrix specified in (|4ip . The quantities oj* and r* can be immediately written in terms of the dimensionless 
variables (l44|) . One can check that both the values of e and |$o| are rescaled out from these expressions. 



APPENDIX E: APPROXIMATE SCALING OF THE SOLUTIONS 



The numerical solutions exhibit an approximate scaling with e, defined after equation (I34p . To see this, we rescale 
the amplitudes and the masses of the flat directions by two different constant factors, 

{$ , i } -> {7*0, 7$o} , {m, m} -f {urn, fim} , (El) 
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We already showed (see for instance the discussion after ea. (fT2")) ) that, under the rescaling (|E1[) , a background solution 
$ (t) , $ (t) is mapped into the solution 7 $ (/i t) , 7 $ (/_t i) . Therefore, /j only affects the timescale governing the 
dynamics of the flat directions. From this, we can see that, under the rescaling (|E1[) . each quantity on the left hand 
side of (f3"4"| is multiplied by /j, while the quantities at the right hand side by 7 . Therefore, under this rescaling, 



(E2) 

To understand how particle production changes under (|E1[) , it is useful to write the corresponding equations in the 
form (|D2|) . Particle production takes place at momenta comparable to the flat direction masses. From the definitions 
(|D3|) , and from eqs. (f3"8")> . we see that 



u u = 0(l) , lo 2M =0(e) (E3) 

For this reason we define 

Wl = W*l , W 2 ,3 = (F4) 



so that all u>i are of order one. We also note that they do not change under the rescalings ([El 
Concerning the matrices I* and J», we have instead (cf. the definitions (|D4j) and eqs. (|4Tj) ) 



h 3 * = -hi* = O (e 1/2 J , 7 23 * = -J32* = O (e) 

Ji3» = J31* = O ( e 1/2 ) , J 23 * = J32* = O (e) (E5) 



while all the other elements are of higher order in e and can be neglected. As for the frequencies, we define 

j _ Il3* f _ I23* j _ Jl3* j _ J23* /Ml 

713 = 7J72 ' 723 = — ' 3 = TV 2 " ' 23 = ~T~ ^ ' 



so that all of these quantities arc of order one, and unaffected by the rescalings QElj) . 

Finally, we rescale the time variable as 77 = 77* /e ; once written this way, the time dependence of the background 
solutions is unaffected by the rescalings (|E1|) (namely, <E> (77) — > 7 $ (77) , $ (77) — > 7 $ (77) ) . This guarantees that the 
three quantities dCdi/dfj are of order one, and unaffected by the rescalings (|Eip . 

In terms of these variable, eqs. (|D2|) read 

= iu)\Pxj - e 1/2 (I13 /3 3 j + J\3 a 3 j) + e 
/ 1 - ~ ~ ^2 

£ a 2 j = ~ e N ^2 Q!2j + ^23 «3j + ^23 P3j ~ $2 J 

e $y = 6 ^ ^2 /?2j - -^23 03j ~ J23 Ot3j + &2jj 

e (3' 3j = e 1 / 2 (J 13 /3y - J i3 a u ) + e f i w 3 ^ + ^23 Ay - J23 a 2j - + a 3j J (E7) 

where j = 1,2,3, and where now prime denotes derivative with respect to fj . As we mentioned, all of the quantities 
Co, I, J, Co' in these equations are of order one, and do not change under the rescalings (IE1[) . In contrast, e is small and 
changes according to (|E2[) . Due to this hierarchy, we can neglect the last term in the first two equations (|E7[) . The 
numerical solutions also show that, as soon as the quanta are produced in an appreciable number, |eay| <§C \u>i ay | 
and |e/3jj <C \u>i /3y|- Therefore, the first two equations of (|E7|) can be written as 

~ -io>i (V 1/2 Oy) "/l3a 3j ;- Ji3/?3j 

~ 7 wi (e" 1 / 2 /3y ) - Z13 /3 3j - J13 a 3j (E8) 



: 1/2 (7i3 ay - J13 /3y ) + 6 ( -i W 3 a 3 j + 723 a 2 j - J23 fcj + ^TT- /?3j 

z w 3 
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while the remaining four equations can be cast in the form 



a 2i = —i w 2 a 2j - h 3 a 3j - J 23 3 j + 02j 

ZUJ2 

I ^2 
02j = i ^2 P2j - h3 P3j - J23 «3j + a2 3 

a' 3 j = I13 (e~ 1/2 aijj - J13 (e~ 1/2 -iQ 3 a 3j + I 23 a 2j - J 23 02j + ^rr- p3j 

0' 3j = lis (V 1/2 0ij) - J13 (V 1/2 ay) + » w 3 foj + I33 P2 3 - J23 a 2j + ^- a 3j (E9) 

From a quick inspection of eqs. (|E8[) and (|E9|) we conjecture that, if we compare the particle production in several 
cases which differ from each other only on the value of e, we will find that 

otij,0ij oc e 1/2 , a 2 j, 02 j, a 3j , 3j oc e° (E10) 

Recalling that the occupation number of the ith mode is, rij = Y) , this implies 

TT-l CX £ , 7T-2, ^3 OC 6° (Ell) 

Equivalently, if we compute the particle production for some given configuration, and we then consider a second 
configuration related to the first one by the rescalings (jElj) . we conjecture that 



{aij , , a 2 j , /3 2 j , a 3j - , 3 j} -> | y ' y ~~ ft? > a2 J > ^ 2 J > a3 J > 3 i 

{m , n 2 , n 3 } -* i "2 , n 3 j (E12) 

The numerical computations that we have performed confirm this scaling behavior (see Section |V| . It is not hard to 
verify that, under the rescalings (|E1[) . the dimensionless quantities defined in (|44|) scale as 

z? i? n n V V ~ \ V 2 \ \ ^ ~\ 

-r* ► -r* , G* — * G* , m* —y—m, — ► — m» , A* — * — ^ A* , A* — » — ^ A» 

— — [background] — > - — - — [background] , fc* — > — fc* , n.j — ► — rij (E13) 
a 77* 7«t), 7 7 

where "background" stands for any background quantity. Therefore, the ratio between the energy density of the 
produced quanta and that of the flat directions, given in eq. (Hoi) , scales as 

Vod -* -9 ^prod (E14) 

7 

Equivalently, we can say that r pro d oc e 2 . 
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